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Linearized perturbations of a Schwarzschild black hole are described, for each angular momen- 
tum £, by the well-studied discrete quasinormal modes (QNMs), and in addition a continuum. 
The latter is characterized by a cut strength </(7>0) for frequencies ui = —17. We show that 
(a) (7(7|0) oc 7, (b) g(r) = at Y = (i+2)!/[6(^-2)!], and (c) q(-y) oscillates with period ~ 1 
(2M = 1). For £ — 2, a pair of QNMs are found beyond the cut on the unphysical sheet very close 
to F, leading to a large dipole in the Green's function near F. For a source near the horizon and a 
distant observer, the continuum contribution relative to that of the QNMs is small. 

PACS numbers: 04.30.-w, 04.70.Bw, 04.20.Jb, 11.30.Pb 
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I. INTRODUCTION 

Gravitational waves propagating on a Schwarzschild 
background probe the nontrivial spacetime around the 
event horizon. If and when detected [1], their signature 
may confirm that black holes exist. Because waves escape 
to infinity and into the horizon, the system is dissipative, 
and described by its spectrum in the lower half frequency 
plane — of interest both for signal interpretation and in its 
own right. The discrete quasinormal modes (QNMs) have 
been thoroughly studied [2-5]; this paper characterizes 
the continuum, about which little is hitherto known. 

For a black hole of mass M (below c = G = 2M = 1) 
and each angular momentum £, the radial functions tp 
of scalar (s=0) or electromagnetic (s=l) waves and ax- 
ial gravitational perturbations (s=2) are governed by 
a generalized Klein-Gordon or so-called Regge- Wheeler 
equation (KGE or RWE) [d 2 x + uj 2 - V(x)]ip(x,u>) = 0; 
x = r + ln(r— 1) is the tortoise coordinate and r the cir- 
cumferential radius. The potential 



F(r!=(l-i 
r 



£(£+!) 1-s 2 



(1.1) 



describing scattering by the background [6], behaves 
as V e (x^oo) = V(x) - £(£+l)/x 2 ~ 2£(M-1) lnir/2; 3 
and V(x->-oo) ~ [t(t+l) + 1 - s 2 ] e^ 1 . We im- 
pose outgoing-wave conditions (OWCs) ip(x— >±oo, uj) ~ 
e iw|x| ^p or x _^ _ 0Oj waves thus go into the horizon.) 

Polar gravitational perturbations are governed by the 
Zerilli equation (ZE) [7], which is the KGE with the po- 
tential V{x) = V(x) + 2d x W(x), where 



W(r) = r + 



3(r-l) 
r 2 (2i/r + 3) 



(1.2) 



F = (£+2)!/[6(^-2)!] = fi/^+l), and v = \{£-l){£+2). 
In fact, W[x) = -g'(x,iT)/g(x,iT) [cf. below (C5)]; 
g{x,uj) is denned in general below. The solutions -0 of 
the RWE and ip of the ZE are related by "intertwining" 
or supersymmetry (SUSY): ip(x,uj) = [d x + W(x))ip(x, uj) 
[2, 8, 9]. Thus, also the two continua are closely related. 

The signal ip(x,t>0) depends on {ip(y, 0), ip(y, 0)} 
through the Green's function G(x,y;t) — f(d<jj/2ir)x 
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FIG. 1: Fourier inversion of G; crosses denote poles in the 
oj-plane. (a) Conventional contour involving a cut on the NIA. 
(b) Modified contour detouring around the unconventional 
pole near F (contour on unphysical sheet shown as broken 
line). For simplicity the mirror contour is not shown. 



G(x, y; u)e~ luJt . Closing the contour in the lower half 
w-plane (Fig. la) separates G into (a) the large semicircle 
(|w| — * 00), giving a prompt signal propagating directly 
from y to x and vanishing after a finite t [10, 11]; (b) the 
QNM poles, giving rise to a ringing signal dominating at 
intermediate t [12, 13]; and (c) our main focus, the cut 
on the negative imaginary axis (NIA) uj = —47: 



AG(x,y; 



-17) = 
-i-y) 



G+(x,y;-ij) - G-(x,y;-ij) , (1.3) 



where G±(— 17) = lim e |o G{— Z7±e) are continuations 
from ±uj > 0. The physical sheet for G+ (G_) lies to 
the right (left) of the NIA. [However, the opposite sides 
are unphysical only for the conventional choice of the cut, 
as in (1.3) and Fig. la.] The continuum is given by AG, 



2 



and for 7 — ► causes the late-i behavior [10, 12]. 
In general, the Green's function G(x, y) — G(y, x) is 



G(x,y;u) = 



V < x , 



(1.4) 



where / (g) solves the KGE with the left (right) OWC, 
and J(g, /; lo) — gf — fg' is their Wronskian. Although 
G is normalization-independent, for definitcness wc de- 
fine f(x^—co,u) <~ 1 • e~ tux and g{x^co,ijj) <~ 1 • e tuJX 
(adopted for Imw > 0, and continued to Imut < 0). At 
a zero of J, / oc g satisfies both OWCs and defines a 
QNM. These are well understood, and we turn to cuts. 

If V has its support in say [— d, d], the OWCs can 
be imposed at ±d, so the KGE is integrated over a fi- 
nite distance; hence, /, g are analytic in ui. The same 
holds if V decays faster than exponentially. If how- 
ever V{x^>— 00) <~ ^2k v k eXkx (here A = 1), typically 
/ has poles ("anomalous points") at ui n = —inX/2 (read- 
ily shown by Born approximation, i.e., a power series 
in e x ). These are removable by scaling f(u>) 1— > x( w ) = 
(w—cj n )f(u>), leaving G unaffected [14]. However, {vk} 
could conspire to "miraculously" make some ui n non- 
singular. For the RWE, n = N = 2r (= 8 for £ = 2) is 
miraculous. Although miracles can be studied by finite- 
order Born approximation, here an exact solution exists 
at T [cf. (1.2) and (C3)], the only miraculous point (for 
any £) [15]. 

On the other hand, for x — ► +00, the centrifugal bar- 
rier does not scatter and has no real effects [10], so we 
should consider the next asymptotic term 



In r r°° 

V e (x^oo) 3-=/ dA(3-2 7 E-lnA)A 2 e- Aa: (1. 

x Jq 



5) 



(7B is Euler's constant). The superposition of exponen- 
tials spreads the poles at —inX/2 into a cut due to the 
power-law tail. In G in (1.4), only g is discontinuous and 
we are led to study Ag(-i^) = g+(— i-y) — g~(— cf. 
(1.3). Section II analyzes this in terms of a position- 
independent cut strength 5(7). The result is checked 
against (a) the 7 ^ limit, and (b) the zero and known 
slope at 7 = T [15]. 

Section III deals with the Green's function, in par- 
ticular AG and the limiting function AG L describing 
propagation from near the horizon to infinity. Of course 
AG oc (7(7) has a zero at 7 = T; surprisingly, there is a 
large contribution, approximately a dipolc, near T, at- 
tributable to a pair of nearby QNMs on the unphysical 
sheet [16]. These poles are absent in q itself, and point to 
a relation to / — in particular to J in (1.4), which must 
have a zero at these positions. This leads to an analytic 
treatment, by linearization about T. Furthermore the 
7^0 behavior is examined to recover G(x, y; t— >oo). 

Some elements of the analysis are computationally 
hard (CH), in a precise sense: consider the evaluation 
of g(x, —i'f) given that g(z— >oo, —17) ~ e 7Z . If the OWC 
is imposed at x = L, one needs accuracy e~ 2lL to ex- 
clude an O(l) admixture of the wrong solution. If V 



does not have finite support (L — > 00) no finite accuracy 
suffices, defeating direct integration of the KGE [17]. In- 
stead, one must continue from the upper half w-planc 
(where growing and decaying solutions are interchanged) 
to the lower one — implicit in all analytic formulas, e.g., 
the Born approximation or Leaver's series (see below). In 
contrast, evaluating a decaying function given its asymp- 
totics, e.g. g(z, i-f) ~ e~ 7Z , is not CH. Calculating g(x, u) 
from {g(z, 10), g'{z, lu)} is also not CH if \x— z\ is finite. In 
Sections II and III, handling the CH parts is the technical 
issue. For z — > — 00, the potential tail is exponential, and 
the calculation is not CH when matching to a finite-order 



Born approximation instead of to f(z, —ij) 



[17]. 



II. CUT STRENGTH 

A. Definition and properties 

On the NIA, g± ~ 1 • e lx satisfy the same RWE. Hence 
Ag <~ • e lx is the small solution oc g(+ij). Since 
g{— uo*) = g*(<jo), Ag is imaginary, so wc introduce the 
real cut strength q [12, 15], to be studied numerically: 



Ag(x, -ry) = iq(j)g{x, +27) . 



(2.1) 



Since g is defined by the OWC at x — > 00, (2.1) de- 
fines q independent of V(x) at any finite x: if, say, 
Vi(x>L) = V2(x>L), the corresponding q\ and qi are 
identical. Leaver [12] has given a formal expression for q, 
which is however nontrivial to evaluate. 

As a result of the SUSY relationship between the RWE 
and ZE, the latter's cut strength q obeys [15] 



(r + 7)g(7) = (T-7)9(7) 



B. Numerical evaluation 



(2.2) 



Since g(x, is not CH, it is simply integrated from 
large x. We compute g(x, —i^±e) by Miller's algorithm 
(failing at e = 0), taking e J, in the difference for Ag; 
see Appendix A [18]. The solid lines in Fig. 2 show the 
ensuing q(j) for £ = 2, 3, 4, and are key results. 

For 1 = 2, 9(7) = at 7 = 0.75, 1.35, 1.90, 2.44, 
2.96, 3.48, 4.00, 4.51, 5.03, 5.54, . . . , suggesting that the 
spacing approaches \. Indeed, by a WKB analysis [19], 



9(7) - 4cos(2tt7) + 0(7~ 1/2 ) , 



(2.3) 



where only the correction depends on £. The zero at T 
and moreover q'(T) w —37.6 agree with [15] 



9'(r) = - 



457T 



137438953472 
37.45 . . . , 



(2100027e 8 + 30148389005) 

(2.4) 



confirming both (2.4) and our numerical accuracy. 
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FIG. 2: Plots of g( 7 ) for (a) t = 2, (b) I = 3, and (c) i = 4. 
The solid (broken) lines are the numerical result (analytic 
approximation) . 



C. Analytic approximation 

Small 7 relates to late-i behavior, controlled by 
(a) many finite-r scatterings, whose effect vanishes ex- 
ponentially, and (b) a few large-r scatterings, which 
therefore dominate. In terms of u{r) = y/l — l/rip(r), 
the RWE reads d 2 u + Pu = [21], where P(r) = 
[r/(r-l)] 2 [o> 2 + r- 3 - §r~ 4 - V(r)}; the ZE follows in 
complete analogy. For large r, both P and P can be ap- 
proximated by P (r) =uj 2 + 2u 2 r~ 1 + [3uj 2 - ^+l)]r~ 2 , 
with errors P(r) - Po(r) <~ b 3 r~ 3 , P(r) - Po(r) ~ b 3 r~ 3 . 
For s=2, b 3 = Au 2 -2v+2 < for small oj, so b 3 = 63+6/V 
has the smaller magnitude. Thus we regard the equation 
with Pq as an approximation for the ZE, and obtain re- 
sults for the RWE through SUSY (2.2). 

Re-expressing the above in terms of ip, the lead- 
ing Vi{x) ~ x~ 3 lna; is reproduced exactly. The error 
is <~ x~ 3 , apparently down by only a factor In a; ~ lni, 
but a pure x~ 3 term does not generate a late-i tail to first 
Born approximation [10]. Thus the leading correction is 
an extra power of x^ 1 <~ t^ 1 (up to logarithms). 

The equation with P is the hydrogen problem with 
fractional angular momentum. The OWC for r — > oo 
selects its confluent-hypergeometric solution [22] and the 
branch of V —u 2 , giving u(r, lo) k, (— 2iujY^ %u} T cr & %U)T x 
U(a-iw, 2a; -2iur) with 2a = 1 + ^{2£+l) 2 - 12w 2 . 

Since we consider 7 — > at fixed r, — iwr w 0, so 



9(7) = 



Au(r, — 17) 
m(r, 17) 



-2 



r(cr+7) sin[7r(o-+7)] 
r(a- 7 ) (2 7 ) 2 ^ 



(2.5) 



For I = 3, 4, (7(7) is similarly oscillatory, but the range 
of validity does not include T (= 20 and 60 respectively) . 



with T the Gamma function (not the special frequency). 
This approximation is shown by broken lines in Fig. 2. 
In particular, it shows that q'(0) = q'(0) = (-1) £ 2tt [10]. 



III. GREEN'S FUNCTION 
A. Evaluation and general properties 

Using J(g-,g + ; —i"f) = —2i-fq(j), one finds 

AG(x, y; -i 7 ) = -2» 7 g(7) ^^TT ' 

with J± = J(g±,f). We evaluate g± as in Section II B 
and use Jaffe's series for / (Appendix B). Fig. 3 shows 
—iAG(x, y; — i 7 ) for 1 = 2; the two contrasting cases re- 
veal a strong x- and y-dependencc. 

The above is the clearest, but more efficient is to 
rewrite (1.4) as (suppressing u — —i-f±e, e^0): 



For I — > 00, one has q'(T) <~ — V96-7T 3 [20], of the same 
order as (2.4). Thus, in each of the analytically accessible 
regimes 7 J. [below (2.5)], 7 — > T, and 7 — > 00 [cf. (2.3)], 
q(j) has a finite limiting behavior for large £, consistent 
with no new features developing in these curves. 



G(x,y) 



\m/m \g(x)/g(z)] 
f'(z)/f(z)-g'(z)/g(z) 



(3.2) 



In principle any z will do. However, if ±z 3> 1, both 
L{z) = f'(z)/f(z) and i?(z) = g'(z)/g(z) approach ±iw, 
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FIG. 3: Plots of -iAG(x,y; -i-y) for £ = 2 and (a) a; = 0.2, 

y = o.i, (b) z = l.o, j/ = -l.o. 



so the exponentially small denominator is hard to evalu- 
ate. Thus choose z <~ 0(1) and (a) compute R(z) (Ap- 
pendix A). This CH step is carried out only at one z. 
(b) Calculate / as above for L(z). (c) Up to irrelevant 
normalization, f(y) and g(x) follow by integrating the 
KGE over finite distances, yielding G. (d) Let e J. 0. We 
have verified that the result agrees with (3.1). 

As said below (1.3), AG(7«0) governs G(t— >oo). Since 
7 4 relates to x — > oo [10], one can neglect Vg for all 
factors in (3.1) except q(j) w (— l) l 2irj (Section II C), 



giving f(x,-ij) 
(-2jx)- e e'< x 



oc x 



e+i 



and g±(x,- 



This yields AG 



(— xy) i+1 -f 2l+2 , elegantly reproducing [10, 12] 

(2£+2)\ i-xyf 



■ii) w [(2£)\/£\] x 
47ri[(2€+l)!!]" 2 x 



G(x,y;t) 



'[(2^+l)!!] 2 t 2t + 3 



(3.3) 



A potential tail Ve(x) ~ x _0! lna; generically leads to 
G <~ t-( 2 ^+ a ) lnt, but in some special cases (including the 
RWE), the In a; factor does not lead to a In t factor [10]. 



B. Limiting x and y 

The physically important limit is —y, x — > oo. Because 
f(y, —27) ~ e 1V and / in general also has an outgoing 



FIG. 4: Plots of -iAG L (-ry) for 1 = 2. (a) For < 7 < 5.6. 
Insets show some regions with two close zeros, (b) The region 
7^4 expanded. The solid line is the numerical result, and 
the broken line is a fit to (3.7) with a 2 = -0.0227. 



part to the right [f(x, —i-y) ~ e JX \, there is a steep posi- 
tion dependence. This is simply a result of the long signal 
propagation time, and is removed if t is measured from 
the first arrival at to(x, y) = x — y. [Closing the contour 
in the upper w-plane readily shows that G(t<t ) = 0.] 
Thus, consider 



G L (t r ) 



lim G(x,y;t +t r ) 

-y—nx 



(3.4) 



with Fourier transform G l (oj) = J(g,f;u>) 1 by (1.4) 
and the normalization of f,g. The cut is 



AG L (i') 



c?7 



AG L (-ij)e 



-it' 



AG L H7) = A[j( 5 ,/;-»7)- 1 ] 



(3.5) 
(3.6) 



[cf. (1.3) for A]. Results are shown in Fig. 4. 

Although AG L (— iT) = 0, surprisingly AG L is largest 
near 7 = T, where it is approximately a dipole (Fig. 4b; 
the dominance is less pronounced if the source is dis- 
placed from the horizon, cf. Fig. 3). The solid line in 
Fig. 5 shows AG L (i') [obtained by integrating (3.5) up 
to 7 — 5.60, hence accurate except for very small t']. 
The cut AG L (— 27) vanishes at (a) the zeros of q(-f) 



[cf. (3.1)], and (b) 7 



3 

2' ' 



The former depend 



■5 




FIG. 5: Solid line (a): continuum contribution AG L (t') for 
1 = 2. Dashed line (b): contribution of the pair of nearby un- 
conventional poles. Dot-dashed line (c): corresponding con- 
tribution from the first QNM, scaled down by 10 3 . 



with the QNM multiplet found analytically to branch 
off from r at infinitesimal a, one concludes that one ad- 
ditional multiplet has to emerge (as a increases) near 
cj = —iT. Rather than the possibilities contemplated in 
Rcf. [15], we speculate that this multiplet may be due 
to the unconventional poles discussed here splitting (as 
they must when spherical symmetry is broken) and mov- 
ing through the NIA as a is tuned [24]. 

Thus it is advantageous to consider another Fourier 
contour going into the unphysical sheet and detouring 
around ui± (Fig. lb), including them as QNM contribu- 
tions (line b in Fig. 5). This slightly reduces the con- 
tinuum (due to the integral along the NIA and often 
neglected as "background" ) . More importantly, suppose 
these poles move through the NIA as a parameter (say, a) 
is tuned, becoming conventional QNMs. In the assign- 
ment of Fig. lb, the total QNM and continuum contri- 
butions are separately continuous; for the conventional 
contour, each is discontinuous. 



only on V(x^+oo) [cf. below (2.1)], the latter only on 
V(x^— oo), scaling with A. If the two tails are sepa- 
rately adjusted the sequences are independent, but for 
the present A = 1, sequence (a) acquires one integer 
member at 7 = T. However, this member simultaneously 
disappears from sequence (b) [cf. above (1-5)], consistent 
with the first-order zero at 7 = T in Fig. 4. 

Some members of the two sequences are close, cf. the 
insets in Fig. 4a. The many zeros also ensure that AG L 
is generically small. (Section IIIE contains quantitative 
estimates.) 



C. Poles on unphysical sheet: numerics 

The behavior of AG L (— 17) for 7 w V can in fact be 
attributed to a pair of nearby QNM poles u>±. Suppose 
G±(lu) « (ai±ia2)/(iw— 4— b±ic), where ai, a 2 , b are real 
and c > so u>± are on the unphysical sheets. In terms 
of £ = 7 — r, on the NIA one has 



AG L H 7 ) 



2ia 2 £ 



(3.7) 



with a\c + 02b = implementing the zero at 7 = T. The 
broken line in Fig. 4b shows this fit for 1 — 2, yielding 



oj ± + iT = =fc - ib w T0.027 + 0.0033i 



(3.8) 



Further, a plot of \J + (£ — i(T+b))\ vs. C > (not shown) 
clearly supports a simple zero at ( m — c [cf. (3.6)]. 

In summary, we have extrapolated Leaver's series to 
the unphysical sheet, revealing nearby poles making the 
largest contribution to the cut. To the best of our knowl- 
edge, this is the first time such QNMs have been found, 
and these obviously have more effect on the dynamics 
than QNMs on the physical sheet at larger |Imu;|. 

In a broader context, consider the Kerr black hole. By 
comparing numerics for moderately small rotation a [23] 



D. Poles on unphysical sheet: analytics 

Interestingly, the above extrapolation can also be car- 
ried out analytically, by assuming the linearization of 
J + {lok— iT) to be valid up to the nearest zero, 



LO + +iT . 



j;(-*r) 



(3-9) 



Since AJ(— iT) vanishes, one does not have to indicate 
the sheet in the numerator. Following the methods and 
notation of [15], one readily obtains 



J(-iT) 



2vN 



2v + 'A V75 
700009 



— - 1 



917504 

with the constants 

N-2 Nj+1 



for I = 2 



73 



-y 

28-11093 _ 8 
3 6 

75 = 



+ 3(2^+3) 



N 



N 



(N-iy. 



-N 



for I = 2 , 
ivN N+1 e- N 



(1 - 2v){v + 1){N - 2)! 
2 24 _ a 



(3.10) 



(3.11) 



(3.12) 



135 



for e 



characterizing f(—iT) and g(—iT) respectively (cf. Ap- 
pendix C). The peculiar integer in the numerator of 
(3.10) has already surfaced in Ref. [15]. 

Calculating J' + {— iT) is formidable (Appendix C), but 
since u) + + iT in (3.8) is mostly real, one can first evaluate 



6 



only Re J' + — much simpler due to the zero in Ag: 
Re J' + (-tT) = |AJ'(-*r) = -y(T)J{g^T),f(-iT)} 

- 4i/(2i/ + 3)' [6A6) 
Ignoring the imaginary part, substitution yields 

= -0.0328 ... for £ = 2 



(3.14) 



comparing favorably with the numerical —0.027, where 
the latter value, also found by extrapolation from the 
physical sheet, need not be more accurate. Taking (3.13) 
as an approximation for J' + provides a bound: given that 
(3.14) has yielded an ui + close to — iT, including lmJ' + 
will render lo + even closer. Indeed, inserting into (3.9) 
the exact (C20) for J' + (-iY,l=2) yields 



UJ + + iT : 



-0.03248 + 0.003436i 



(3.15) 



as our final estimate; especially the agreement of Imw + 
with (3.8) is remarkable. 

Combining (3.5) and (3.7), our linearization of 
J immediately yields ai « Re( J' + {— ir)) _1 sa 
-Reu + /J(-iT) ss -0.0426. The agreement with the 
numerical fit 012 ~ —0.0227 is slightly worse than the one 
for lo + itself. 

For £ = 3, the same procedure would lead to lo + + iT w 
—0.847 + 0.0433i — outside the radius of convergence of 
the series for J + (uj), which is | due to its poles at the 
anomalous points of /. Heuristically removing these 
poles by rescaling J, e.g., by linearizing H(u>) = J(ui) x 
sm(2wioj) / (ui+iT) [cf. (B4) below], does not change the 
prediction for uj±. Besides, the fact remains that J(uj) 
has structure on a scale of O(l). Thus, already for £ = 3, 
currently there is only slight evidence for ui± near — iT. 
Similarly, for large £, the asymptotics of q'(T) given at 

the end of Section II B yield (3.14) as uj + + iT 0.051I\ 

far outside the permitted range of extrapolation. While 
we have not fully studied the large-£ asymptotics of 
Im J' + (—iT), this is unlikely to change the conclusion that 
our calculation is valid for £ — 2 and possibly £ = 3 only. 



E. QNMs versus background 

In terms of t', also the QNM contributions can be ex- 
pressed for x, — y — > 00 in a position-independent man- 
ner. The various magnitudes and phases arc then com- 
pletely determined — unlike the usual situation in (Q)NM 
analysis. Fig. 5, line c shows the contribution of the 
first QNM (i.e., the pair with the smallest damping, at 
u>i = ±0.747- 0.178i), reduced by 10 3 . Insofar as the 
low QNMs arc likely to be the most important for signal 
analysis, the comparison in Fig. 5 provides the estimate 
that the "background" is small. 



The above numerics refer to the most important case 
£ = 2. The regions near the much larger T(£>3) are 
currently unattainable numerically. In any event, Sec- 
tion HID shows that the unconventional poles, if present 
at all, must be further removed from T, so their influence 
on the cut is likely to be less pronounced. 



F. Zerilli equation 

The corresponding quantity for the ZE is trivially 
AG L (— 27) = AG L H 7 ), by (3.6) and J — J [25]. Thus, 
the discussion of Fig. 5 also applies to the ZE. 

Since the ZE has a QNM at -iT but the RWE docs 
not [15], it might appear that the two G L 's should have 

different pole structures there. Indeed, G does have a 
pole, with residue oc f(y, —iT)g(x, —iT). However, — iT 
is anomalous for the ZE (but miraculous for the RWE), 
so / does not contain the normal growing solution, elim- 
inating the residue in G L . 

Physically, for a source not at the horizon, the axial 
(RWE) dynamics does not have a term e~ rt , but the po- 
lar (ZE) dynamics does — with, however, vanishing am- 
plitude as the source approaches the horizon (y — > —00). 



IV. DISCUSSION 

Like the hydrogen atom in quantum mechanics, the 
Schwarzschild black hole is the simplest compact object 
in relativity, and its spectrum also contains a continuum 
in addition to the discrete part. We have characterized 
this continuum, recovering the behavior both for 7 = 
iu — > and near the miraculous point T. The present 
numerical method cannot access much larger 7, but the 
pattern already appears clear. 

Although AG L (— iT) = 0, it surprisingly is largest 
near T. Moreover, for a limiting source and observer, 
the cut contribution is small relative to the QNMs — of 
relevance once gravitational waves are detected. 

Signal analysis in terms of a V(x) is formally an inverse 
problem. These are well-studied for closed systems [26]: 
two real spectra determine V on a finite interval. For 
open systems, we conjecture that one complex spectrum 
suffices — provided the discrete QNMs are complete. The 
cut (rendering them incomplete) is likely to hamper in- 
version, but, intriguingly, the extended family of QNMs 
(including the unconventional ones) may be complete for 
the dynamics (apart from the prompt signal) and permit 
inversion, even though V is not finitely supported [13]. 
The present at least shows that one pair of nearby poles 
on the unphysical sheet already dominates the cut. 

These questions may be explored through solvable 
models with potential tails [3, 27]. Some aspects of the 
RWE and ZE can also be analyzed asymptotically [19]. 
Numerical algorithms (e.g., generalizing the continued- 
fraction method [12]) valid on the NIA and even into the 
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unphysical region would also be useful, allowing QNMs 
there to be studied directly rather than through extrap- 
olation. 

Finally, it would be instructive to evaluate the con- 
tinuum for a Kerr hole, to study any unconventional 
poles, in particular their movement and possible emer- 
gence onto the physical sheet as the rotation is increased 
from zero [28]. 
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APPENDIX A: LEAVER'S SERIES 

Suppressing the parameter ui, we first calculate an 
wnnormalized g. Define R = g'/g = g'/g. Choose z 
and compute R~[z) by Miller's algorithm of downward re- 
cursion on Leaver's series [29, 30]. With, say, g(z) = 1, 
g(x) is then integrated trivially for all x. Thus, the norm 
TV = lim 

x^oo e %u)x j ' g{x) is not CH; convergence provides 
a check. Finally, g — Afg. This procedure requires only 
one CH calculation [for R(z)], instead of one for each x. 

In principle, the result should be z-independent; in 
practice we take \z\ ~ 1. However, we have checked 
that the R(z) thus obtained obey R' = V — uj 2 — R 2 . 
Such tests show that Miller's algorithm works well at 
least for -4.25 < Imw < -1.25, 1(T 4 < Rcui < 0.1, and 
—0.5 < z < 1.5; by optimizing z, the range on the NIA 
can be slightly extended to 0.5 < 7 < 6.0, varying some- 
what with I. One also verifies that i?(z^±oo) — > ±iu. 

In practice we calculate (7(7) using a small e ^ and 
some definite x in (—5,0); the smaller e, the wider the 
range of permissible x. Convergence for e — > is rapid, 
and the x-independence in particular gives an accuracy 
estimate. On this basis, the decimals given in the main 
text should be significant, and the error bars in the figures 
should be tiny on the given scale. 



APPENDIX B: JAFFE'S SERIES 

This appendix presents Jaffe's series [29, 31], with par- 
ticular attention to the anomalous and miraculous points. 
The series, first used for the H^ ion and applied to the 
present problem by Leaver [29], expresses the outgoing 
solution, /(r|l,w) - e~™ x = {r-l)~ iu - 'e'^ ', as 



The coefficients a n {uj) satisfy 

a n a n +i + P n a n + J n a n -i = , (B2) 
n = 1, 2, . . ., with a n< o — and (s — 2 throughout) 
a n = (n+l)(n+l— 2iu>) , 

(3 n = -2n 2 + (8iuj-2)n + 8u 2 + Aiu - £(£+1) + 3 , 
7„ = n 2 - Aium - Au 2 - 4 . (B3) 

The normalization ao can be chosen freely (see below). 

Jaffe's series leads to some important insights. First, 
a n has a zero at oj = — i(n+l)/2. So at to n = —in/2 
(n = 1, 2, . . .), a m > n and therefore f(u>) generically have 
simple poles, removable through scaling by uj — ui n , in 
effect killing the a m<n and making a m >„ finite [15]. Then 
the leading behavior is x( r > w n) ~ (r—l) n ^ 2 ~ e lLJnX , an 
incoming wave. Thus, at these anomalous points, / cx \ 
is outgoing (being the analytic continuation of outgoing 
waves defined for Ivaui > 0) and incoming. 

However, according to (B2), a n (hence, all a m >„) can 
be non-singular even when a„_i = 0, if /3„_ia„_i + 
7n-ian-2 = 0. This is called a miracle; in the present 
case for I = 2, it occurs at n = 8, i.e., at u = —iT [15]. 

Because of the poles at the anomalous but not at the 
miraculous ui n , we choose to calculate \ with 



a a = 



sin 2777 



(B4) 



in which case all functions are nonsingular at any u n 
(both anomalous and miraculous). However, if needed 
[e.g. in (3.6)] we can always find the normalization by 
taking r | 1; cf. Appendix A. 

Finally, note that f(u>) has no cut in the w-plane, as is 
also evident from more general considerations. 



APPENDIX C: FREQUENCY DERIVATIVE OF 
THE WRONSKIAN 



1. Preliminaries 

Writing f± = d u f(u))\_ iT and similar for g, we are 
interested in (arguments — iT suppressed where possible) 

J'+ = 9i+d x f + gd x fi - f\d x g - fd x g 1+ . (CI) 

The differentiated RWE 

[dl-T 2 -V{x)\h{x) = 2iYf{x) (C2) 

etc., plus nontrivial differentiated OWCs, determines f\ 
and gi + . The calculation proceeds in r, using x(r) as a 
shorthand only. One introduces the SUSY generator 



/(r,w) = (r-l)" 



-itj Jliuj Aujt 



n=0 



r-l 



(Bl) 



(C3) 
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and the secondary solution £2(r) = £i(r) J^dt [t/(t—l)] x 
S,i 2 (t), where the integral is elementary [15]: 



aw 



i 

75 

( 

+ - 



,N(r-l) 



N-2 

E 



[N(l- 



J 



1 



3=0 

Nr (r-l) N - l [2vr 2 - {2v+'i)r + 6] 



Av 2 N{2vr + 3) 



(C4) 



Note that the decay of £ 2 near the horizon is due to a 
high-order cancellation, cf. the first line of (C4). This 
behavior complicates the entire calculation, right down 
to the floating-point evaluation of (C20) or its higher-^" 
counterparts. Namely, for £ = 2 (£ = 3) the first 3 (16) 
digits cancel in the rational and transcendental [largely 
oc Ei(iV)] contributions to Im J' + , so that very high accu- 
racy is needed already for moderate I. 
The Wronskian is readily evaluated as 



1 . 



(C5) 



With £j = £i + 7j?2, the outgoing functions are / = 
£3/(21/ + 3) and g — 2^iV£s/75, while the incoming ones 
read f(iT) = (2z/ + 3)7V£ 2 and g(iT) = fi/(2i/); the latter 
has been used already in deriving the final line of (3.13). 



2. Differentiated wave functions 

Using a Green's function approach, (C2) is solved as 
iN 



fi(r) 



2i/ + 3 



f 



d(r)/ dt 



& (*)&(*) 
t-i 

*&(«)£»(*) 



t-i 



(C6) 



where the undetermined integration constants reflect the 
possibility of adding a homogeneous solution to (C2). Of 
course, /1 itself is not arbitrary, and the second of these 
constants follows from fi(x— >— 00) ~ [— ix + 0(e x )]e~ Tx . 
The required asymptotics of the integrals in (C6) are 
straightforward; for definiteness we write the answer as 



^iV- 1 (2 ^ /+3)/l = + N^x + 73D - N~ 2 } 
-&{C + l3 {B+N- 1 x}+p}; 

t 



B(r) 
C(r) 



dt- 



f dt 



t- 1 

t-i 



Mm 2 (t)-N- v 
D{r) = [ 



dt 



t-i 



(C7) 



(C8) 



The integrals (C8) can be evaluated by partial-fraction 
expansion, yielding only elementary functions plus 



Ei(r) = J r ^ co dte t /t [32]. However, the full primitives are 
cumbersome — not surprising since already £ 2 is involved. 
Hence, we only give results when required. 

One cannnot find (3 from low-order asymptotics, as 
was the case for 73 in the calculation of / itself [15]. For 
£ = 2, we have obtained the Born series for f(r~l,w) 
by computer algebra. Only after is taken in each 
term does one set u = —iT, comparing the result to the 
analogous expansion of (C7). The first eight terms agree, 
as they must. The ninth [0((r— l) 4 )] terms agree if 

/3 = 9Ei(-8 ) + |( 7E+ ln8) + l^e-«, (C9) 

and for this value the two series coincide. The transcen- 
dental Ei, 7e, and In 8 only occur to cancel their counter- 
parts in the expansion of C (r) , defined with a lower limit 
r = 00 to ensure convergence — the Born series for f\ es- 
sentially involves rationals only. Indeed, Ei(— 8) does not 
occur in (C20) below. 

For general £, we calculate (3 by requiring that 
d 2 [(r — 1) IUJ f{r,w)]_ iT be single-valued near r = 1. This 
generalizes the determination of 73 and hence / in [15], 
by demanding that fi(r)+if(r) ln(r— 1) be a power series 
near r = 1. Using the latter, presently we find that 

f 2 (r) - 2 i N- 1 f 1 (r) + 2tf 1 (r)ln(r-l) 

+ [2N- 1 - ln(r-l)]/(r) ln(r-l) (CIO) 

should not involve ln(r— 1). The equation for / 2 analo- 
gous to (C2) yields [cf. (C6)] 



/ 2 (r)-|/i(r) = 2iJV 



£2 



-Zi{r)J 



t-i 



■ (Cll) 



Insertion into (CIO) using (C7) for /1 throughout even- 
tually yields that all logarithms indeed cancel if 



/3 = 9Ei(-AT) - 



(18i/+9)7s , 273 
l/<-'.\- 



N-l 



N-2 



E a e N E fcl 



N 



75 



kl 2v 2 N 



7E +lniV+^ 
+ Np . (C12) 



3=1 3 L k=j 
The residue p = Tes[r^f(r)B(r)/(r—l)] r —i reads 

N-2 



p = {2v+2,fc N + (9-4z/ 2 ) CjV -i + 9 



v ■ 



(C13) 
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2N 

NP 



- 1 



{2u+i) 2 N 
p!(JV-l)! 

JV-l 



JV-2 



N-i- 1 y-v- 9A^ fe - J 

i(p-i)' 'hh'Jip-iW 

N-l 

E 



(2^+3) 2 iV w - 1 



1 



j=l k=0 

(j_l)!2P+i 
(P+J')i 

JV-l 



AT -2 



For ^ = 2, (C12)-(C14) are verified to reduce to (C9). 



k=j ^' 

(2v+3) 2 N N 
(N-j)\ 



(N+p)\ 

6(2^+3) 



(2 W+P -1) 

N N-1 



(N 



N-l ~\ 



(C14) 



The calculation of gi+ may seem more difficult concep- 
tually, involving analytic continuation in r. It actually is 
easier technically, since the exact form does follow from 
low-order asymptotics. The counterpart to (C6) reads 

^T-Ur)j dt — -a(r)y ' 

(C15) 

where the lower limit for the second term ensures 
gi(r— >— oo) — > 0. We select <?i+ by continuing r from 
— oo to the physical r > 1 through the upper half plane. 
The remaining integration constant is made explicit as 



2ivN 2 



75 x 
N 



+ a 



-€i(r) 



eft 



t-1 



(C16) 



Comparison to the differentiated asymptotics of g(u>), 
gi(r) 3x 3x 3 



3x 
2vr 



0(r- 2 ) , 



(C17) 



ie Nx/2 2vr Nr 2 vNr 

yields a = j 5 /(4n 2 N 2 ) in 0(r°) of the last factor in 



(C17); agreement of 0(r 1 ) provides a check. 



3. Evaluation of the differentiated Wronskian 

Substituting (C7) and (C17) into (CI), the prefactors 
of the t- integrals are all handled using (C5), yielding 



(2^+3)75 
2ivN 2 



75-73 

N 



-&(*)&(*) 



73-75 
Av 2 N 2 



- l5 N- 2 +p + C(r) + ( 73 +75)B(r) 

+ 7375-D(r) + 2 73 iV- 1 x . (C18) 

Differentiating with respect to r, the r-independence of 
(Cf 8) is readily verified. Indeed, the rhs is merely a regu- 
larized j 1 °° dt £,3£,5t/(t— I), which unfortunately diverges 
at both ends. Choosing r — > —00, the first term vanishes. 

Inspecting (C8), it is gratifying that the individual ex- 
ponential divergences (for r — ► —00) in B, C, and D 
indeed cancel in (C18). Taking the i-contours in the up- 
per half plane fixes the contributions of the singularities 
at t = 0, 1, yielding Re J' + as in (3.13). In full, one has 



lim \C(r) + (73+ 75 )S(r) + 7375£>(r) + 2 73 iV- 1 x] 

r^ — oo+iry 



f (5 -2)(^+ln^)-9Ei(-^), 9Ar275 



4 7| ^ :/Ar , , „_,73 / -, 73 



Ei(iV+ir?) + 2tt^ 1 - — 

n V 75 



JV-l 



75 



-A* 



2^ 2 iV 



^ fc! 



fc=j 



373(2^-1) 2(2^+3)(2t/+5)(8z/ 2 +8^-3)7 3 



2 73 (9-6^)73 iV-3 973756^—' < ! 



+ TV 2 + 2^ 2 iV-2 



(2z/+3)iV 2 

AT-4 

J! 



E 



4^ 4 ^ A"+ 3 



9N 2 (N-l)(N~2) 
(12i/+9) 73 iV jV e- JV _ I 



3(2^+3)75 , 2 73 ^ 2 I 973e-_^ y^^ W 



4^2 



fri 3 75 3 tt ifo kk\(j-k)\ v 2 ' 



75 frinKN-jy. 
| 9 73 ^ 2 g^- fc - 2 fc! 



=3 fc=0 



+ 



4^7 3 e" 



-N 



N-4 



75 



973 e 



-JV 



AT-5 



J! 



fc! 



75 



E^ + ^ + 3E^ ; 



i=i 



fc=0 



(C19) 
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with Gj = J2k =] 2 [k\(N-2+j-k)\] \ The infinitesimal 
f] > selects the upper branch of F,i(N+irj). 

Finally, Eqs. (3.11), (3.12), (C12)-(C14), and (C18)- 
(C19) combined yield a tedious, but explicit exact result 
for J' + ; some cancellation occurs. For £ = 2, one has 



J' + = 49 ^38 [-17122265640585(7E+ln8-m) 

- 245810518235861775 Ei(8+^)e" 8 
+ 36326230655979688] . (C20) 
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